Glassy behavior of light 
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We study the nonlinear dynamics of a multi-mode random laser using the methods of statistical 
physics of disordered systems. A replica-symmetry breaking phase transition is predicted as a 
function of the pump intensity. We thus show that light propagating in a random non-linear medium 
displays glassy behavior, i.e. the photon gas has a multitude of metastable states and a non vanishing 
complexity, corresponding to mode-locking processes in random lasers. The present work reveals 
the existence of new physical phenomena, and demonstrates how nonlinear optics and random lasers 
can be a benchmark for the modern theory of complex systems and glasses. 

PACS numbers: 64.70.Pf, 75.10.Nr, 42.55. Zz 



The first marriage between statistical mechanics and 
lasers dates back to their early development 0, 
Since the seventies, many authors have outlined that the 
threshold for lasing can be interpreted as a thermody- 
namic phase transition; and these ideas spread out in the 
field of photonics, later embracing also nonlinear optics 
(for a review see e.g. Q)- In recent articles, "i*, Y, 
the statistical properties of laser light in homogeneous 
cavities has been studied taking into account nonlinear 
phenomena, like gain saturation and intensity dependent 
refractive index. This nonlinearity gives rise to an inter- 
action among the oscillation modes, which in turn pro- 
duces new interesting effects. Specifically, by mapping 
the dynamics in an ordered Hamiltonian problem, the 
authors of Ref. Q predicted a critical behavior (a phase 
transition) of the laser mode-locking process. 

In the latter example, the statistical mechanics of or- 
dered systems is applied to study light propagation in 
amplifying homogeneous non-linear materials. In recent 
years, the locally inhomogeneous character of matter, 
and in particular the disordered nature of these inho- 
mogeneities, are becoming more and more important. 
Specifically, relevant attention has been dedicated to 
light amplification in random media and random lasers 
fg, 10, 11, 12, 13, 14]. This is a fascinating topic that 
bridges various fields like light localization and diffusion, 
thermodynamics, nonlinear physics and quantum optics, 
and it has relevant fundamental and applicative perspec- 
tives, as in bio-physics 15]. Methods of statistical me- 
chanics have not yet been applied to the propagation of 
light in non-linear disordered active media, but the rich 
behavior observed in glasses (aging, memory effects, etc.) 
can be foreseen 0, . 

In this Letter, we analytically study the statistical 
properties of the modes of a random optical cavity. Var- 
ious physical settings are embraced by this problem: for 
example a micro-structured cavity filled by an active soft- 
material, like doped liquid crystals, cavity-less random 
lasers, or even a standard laser system with a disordered 
amplifying medium. In all these cases, the disorder varies 
on time scales much longer than the optical fields (it is 



"quenched" O) and, for any realization of the system, 
the supported electromagnetic modes interact because of 
the nonlinearity of the resonant medium. 

We first show that light propagation is described by 
an Hamiltonian with the pumping rate acting as inverse 
temperature. The mode amplitudes are taken as slowly 
varying and their phases play the role of "spins" in an 
Hamiltonian, which turns out to be a generalization of 
the XY model ^3] and of the so called fc-trigonometric 
model The thermodynamics of the model is stud- 

ied using the replica trick |16|], and a one step replica 
symmetry breaking transition is found pol ]. Our results 
show that when the average energy into each mode in- 
creases (i.e. the "temperature" is decreased), the sys- 
tem undergoes a glass transition, meaning that its dy- 
namics slows down and an exponentially large number 
of states appears, with a non vanishing "complexity" 
|20j ] . They correspond to "mode-locked" states of a ran- 
dom laser. This treatment somehow generalizes to disor- 
dered systems early works on multi-mode lasers (see e.g. 
'2T, '25', "2^ ) , and can be also applied to other problems 
involving multi-mode interactions in a nonlinear medium. 
By providing an analytically treatable statistical model, 
not only we unveil the complex behavior (in the meaning 
of modern glassy physics) of light in active disordered 
media, but we also demonstrate the possibility of using 
them for testing the replica symmetry breaking transi- 
tions. When studying the glass transition of atomic or 
molecular systems, due to the predicted kinetic arrest, 
the interesting time scale becomes extremely long with 
respect to that experimentally accessible. Due to the 
intrinsically fast photon dynamics, it is reasonable to ex- 
pect that the system can be equilibrated much closer to 
the transition, and this should provide new and inter- 
esting experimental data on kinctically arrested "photon 
glasses" . 

We consider a dielectric resonator described by a re- 
fractive index profile n(r). The time dependence of this 
quantity is of no interest here, taking place on a time 
scale much slower than that of the photon propagation. 
Our approach follows the standard coupled mode-theory 
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in the time domain [2J| . The Maxwell's equations in the 
presence of a nonlinear polarization Pnl are; 

V X H(r, t) = eon^{r)dtnr, t) + dtV^dv, t, E, H) 
VxE(r,i) =-/xoatH(r,t) . 

(1) 

The fields E(r, t) and H(r, t) are conveniently decom- 
posed in normal modes corresponding to the solutions 
of the linear problem, Pnl — 0, (frequencies ujm, eigen- 
vectors E,„(r), H„i(r), and "amplitudes" am)- For later 
convenience, they are cast in the form 



E(r, i) = 3?[X;„amE„j(r)exp(-«cj„ii)] , 

H(r, i) = y/uJ:;^ 3fJE„amH„(r)exp(-iw„t)] , 



(2) 



where E„i(r), H,„(r) solve the eigenvalue problem 

^2 \f E„ 

and are normalized such that 
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i ^{£o«'(r)|E„,(r)|2 + ^o|H,„(r)|2}dV^ = 1 



(3) 



where V is the standard normalization volume, over 
which periodical boundary conditions are enforced. With 
these notations, the total energy stored in the cavity is 

In the presence of nonlinearity, the amplitudes 
time-dependent and their evolution is described by cou- 
pled equations that can be derived using standard per- 
turbation techniques [23| and take the general form 



dam{t) 
di 



(4) 



The integral in Eq. Q is extended to the cavity vol- 
ume where the nonlinear polarization is different from 
zero. The quantity Pm(r) is the amplitude of the com- 
ponent of the nonlinear polarization oscillating at uJm, 
'Pnl = ^[J2m'^'mG:>q>{—iLL!mt)] and its explicit expres- 
sion as function of the fields depends on the considered 
non-linearity. For isotropic media the leading cubic terms 
are written, in our notation, as |25| 

-Pm(r) =^Kapjs{t-^7n; OJp,UJq, -UJr, r)E^ E'^ Ef.apaqa*. 

Here Kap^s — ^^p^q'jJrXap-^& and x is the third or- 
der response susceptibility tensor (explicit expressions 
are known for example in the two level approximation 
H^HEHlDi the sum over Cartesian indeces is imphcit. 
The coupled mode theory equations in a nonlinear cavity 
read hence as 



where the sum is extended over all the modes and 



(5) 



■jmpqr ' 
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■ Kap^si^^m; t^p, t^g, -UJr , r)E^^E^ E^ Ef.dV . 
Jv 

(6) 



Under standard approximations |22, |25| , the tensor g is a 
quantity symmetric with respect to the exchange of any 
couple of indeces and can be taken as real-valued, ne- 
glecting processes like self and cross-phase modulations 
between modes (S-XPM). These mechanisms can be in- 
cluded in the model (see below), but will be here ne- 
glected for the sake of simplicity. Introducing the (real 
valued) function 



~7 ^PQ''^9spqr^s^p^q^r 



(7) 



and taking into account radiation losses and material ab- 
sorption processes, represented by the coefficients a„i and 
light amplification through 7™ (see also 01 )i the result- 
ing equation of motion for am{t) is written as: 



da,i 



dt 



dH 

da* 



da*m 

(8) 

where V. = I]m("™ - lm)\am? + H, and having in- 
troduced as usual the noise term (see e.g. 0;|23); with 
{Vp{t)'nq{t')) = 2fcBT'f,at/i(5p,(5(i-t'), weighted by an effect- 
ing temperature T^atin with ks the Boltzmann constant. 
The previous equation is a standard Langevin model for 
a system of N particles moving in 2N dimensions (rep- 
resented by {5Ram, 9am}m=i..Ar) 1^ and its invariant 
measure is given by exp(— H/fcsTbath). Note that the 
latter is not affected by S-XPM terms, as we will discuss 
in future publications and originally addressed in Q. 

We consider a large number of modes in a small fre- 
quency interval ojm^^o pumped and put into oscillations. 
We can write am{t) = Am{t) exp[i(^m(t)], which is useful 
as we expect Am to be slowly varying with respect to 
LPmi as discussed below. In previous works, with refer- 
ence to standard multi-mode lasers, |2ll l2^ the phase- 
dependent terms in © , were always averaged out by as- 
suming the Lpm rapidly varying, independent and uni- 
formly distributed. The resulting equations determine 
the oscillation energy £m into each mode. However, for 
an increasing number of modes, nonlinear beatings in- 
duce non-trivial light dynamics which is mainly due to 
the rapidly varying phases, while the amplitudes can still 
be taken as slowly varying I 211 0, . 

Summing up, the Hamiltonian depending on the rel- 
evant dynamic variables, the phases ipm, is written as 



7Y(G, ip) ^l-lo + T.spqrGspqr COa{ips + Ifip - Ifiq - (fir) (9) 

where Ho = J2m('^m ~ lrn)A^ is an irrelevant constant 

we will con- 
sider these G coefficients as "quenched" (due to the slow 
t dependence of Am), and the relevant phase space is 
reduced to that spanned by Lpm- 

If the cavity is realized by a random medium, as de- 
scribed above, the coupling coefficients G are random 
variables. Their values depend on the mode profiles, the 
resonant frequencies, and on the quenched values of the 
energies Em = "^m^m i'^ each mode, which vary with 
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each realization of the cavity for a given pumping rate. 
For these reasons we take Ggpgr as random Gaussian vari- 
ables. Due to the fact that the mode frequencies are typi- 
cally symmetrically distributed with respect to the tran- 
sition frequency of the amplifying atomic system (and 
correspondingly the signs of the nonlinear susceptibility 
terms largely varies, see e.g. H^H^) and their spread is 
small, LUm'^i^o, G is taken with zero mean value (G) = 0. 
The latter hypothesis can be removed by a suitable, but 
not trivial, generalization of the treatment below [l6l |. 

The Gs roughly scale as {A^)'^V~^^'^ ^ as one can derive 
from Eq. ® by using that is 0(V"-i/2) and K^fj-ys 
is a random variable integrated over V . Recalling that 
uom ^ 1^0, t^o(^^) measures the average energy per mode. 
By a simple rescaling, the invariant measure can be writ- 
ten as exp[-PH{J,(p)], where Jspqr = Gspqr/{go{A'^)'^) 
has standard deviation oc oc l/N^^^ and go is 

a material-dependent constant. Note that this scaling 
of the Js guarantees that the Hamiltonian is extensive 
[T^ |20| . The parameter that controls the phase transi- 
tion is then /3 = l/T = {A^ygo/kBTbath- Thus, in what 
follows transitions obtained as /3 is increased (i.e. the 
adimensional effective temperature T is decreased), can 
be controlled by increasing the amount of energy stored 
on average into each mode (i.e. the pumping rate). 

The calculation of the thermodynamics goes through 
the evaluation of the partition function 



Z{J) = Jd(pe 



(10) 



and, more specifically, of the free energy /(/3) averaged 
over the quenched disorder, that can be written, using 
the replica trick , as 



-/3/(/3) 



lim — logZ(J)= lim lim 



Z"(J) - 1 



By a standard saddle- point e valuation, the replicated 
partition function is Z"{J) = exp[— iVn min ^(q)], 
where the overlap matrix q is defined by qab = 
J2j ^'Hd one has to find the overlap val- 

ues q that satisfy the saddle-point equation. The free 
energy is then given by 



/?/ = - lim(niV)^i{e 

n— >0 



1} 



lim fl(g) = /3(j>{q) , 

n— >0 



where 
Z{q) = 



4 



3/32 



'^^a>b\qab\'^ 



log Z{q) , 



diys" exp 



(11) 



This function is very similar to the one that describes the 
Ising p-spin glass for p = 4, see e.g. 0. The derivative 
with respect to qab of f{q) gives the saddle point equation 




FIG. 1: The Irsb overlap g as a function of the reduced tem- 
perature. The full line is the stable part of the curve, the 
dashed line is the metastable part. Inset: the complexity 
E(r) as a function of the temperature. It jumps discontinu- 
ously at Td where the metastable solution first appears and 
vanishes at Tc where the glass transition takes place. 



where the average is on the measure that defines Z[q). 
In the following, we will restrict ourselves to consider real 
values of qab, without taking explicitly into account the 
trivial rotational symmetry [23. 

The replica symmetric (rs) solution corresponds to 
qab = q and as usual it turns out that g = is the cor- 
rect solution. The RS free energy is simply fas — — /3/4 
as in the Ising p-spin glass. In this regime the are 
uniformly distributed in [0, 27r), inde pen dent and rapidly 
evolving, as originally considered in |2ll l22 | . 

In the one step replica symmetry breaking (Irsb) 
ansatz one divides the matrix qab in n/m blocks of side 
TO 16]. The elements in the off-diagonal blocks are set to 
while in the diagonal blocks RS is assumed and qai, = q. 
In this case, Eq. Hll|l becomes: 

mm,q)= - (/3V4)[l + 3(l-TO)g4_ V] 



TO 



-ilog/o"dzze^*/o(/?Zz)™ , (13) 



where /„ indica tes the modified Bessel function of first 
kind and / = \/2(p. The self-consistency equation that 
defines q is 



J^dz: 



J^dzze-^I^iPlz) 



(14) 



qab = {eM^iv" - v'')]) 



(12) 



These expressions are identical to the Irsb free energy 
for the {p = 4) p-spin model, the only difference being the 
presence of the Bessel function instead of the hyperbolic 
cosine in the last integral. Starting from this expression 
one can repeat the analysis of to derive the full phase 
space structure of the model at the Irsb level. 

In Fig. 1 the solution q of Eq. H14|l for m = 1 is re- 
ported as a function of T. At high temperature q — 
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and the RS solution is recovered. On lowering the tem- 
perature, a solution q ^ first appears at (dashed 
line). However, it becomes stable only below the thermo- 
dynamic glass transition temperature Tc < Td (full line), 
as in standard first-order transitions. For Tc < T < Td 
the phase space is disconnected in an exponential num- 
ber 7V(r) = expiVE(r) of states. The Irsb complexity 
J:{T) is reported in the inset of Fig. 1. At T = Td a 
dynamical transition is expected to take place [l7|, |20| . 

The stability of the iRSB solution can be studied 
(within the assumption that Qab is real). It turns out 
that the Irsb solution is thermodynamically stable for 
all temperatures, so in this case no Gardner transition is 
present [MEllSl- 

The presence of a dynamical phase transition at a given 
value of the random laser pump intensity implies differ- 
ent interesting physical phenomena. These could be ex- 
perimentally investigated by studying (for example via 
heterodyne experiments) the self correlation function of 
a specific frequency (wm) component of the electric field 
in the cavity: 

C{t,UJrn) = ^TuA^m (exp{i[.^™(t + t) - ^rair)]}) ^ ■ (15) 

On approaching Td from above (i.e. on increasing the 
pump power) , the dynamics of phase variable (t) be- 
comes slower and slower and C{t,ujm) is expected to de- 
cay towards zero in longer and longer times. At the 
dynamical transition, the dynamics of the ip^s becomes 
non-ergodic, they are no longer able to explore the whole 
phase space and the function C(t,ujm) decays towards 
a plateau. In other words, the mode's phases ipmit) - 
beside from small oscillation- are locked to some "equi- 
librium" values ("random mode locking"). Due to the 
non-vanishing value of the complexity S at Td, however, 
an exponentially large number of such equilibrium posi- 



tions exist, so giving rise to many different time structure 
of the electric field in random lasers. On further increas- 
ing the pump power, the complexity S decreases and the 
ideal glassy state is reached at T=Tc. Below this value 
the number of equilibrium states is not exponential in 
(E = 0). The equilibrium states below Td are diffi- 
cult to reach because the needed time scale diverges ^3 ■ 
interesting phenomena as aging, memory effects, and his- 
tory dependent responses are expected to take place for 
T<Td. 

All these non equilibrium phenomena arc theoretic ally 
predicted and, to some extent experimentally verified |l7| 
in material systems (structural glasses, spin glasses,..). 
The absence of conclusive experiments is mainly due to 
the long time needed to "equilibrate" real glasses below 
Td, a time which is dictated by the elementary time of 
the atomic/molecular dynamics in condensed matter (ps 
time scale) . Photon dynamics in cavities is at least three 
order of magnitude faster. This observation drives us 
to propose the random lasers in disordered media as a 
benchmark to test experimentally the outcome of those 
theories, as the replica symmetry breaking. 

In conclusion, the multi-mode dynamics in a random 
laser cavity has been investigated by statistical physics 
techniques, and a one step replica symmetry breaking 
phase transition has been found. Our results emphasize 
two important points: i) the light propagation in non- 
linear disordered media shows the same complex behavior 
of the dynamics of glassy systems (aging, memory, ...); 
ii) due to the faster photon dynamics with respect to the 
atomic one, it is possible to use the random lasers as 
systems for experimentally testing the replica symmetry 
breaking transitions. 

We are pleased to thank T. Castellani, L. Leuzzi, and 
F. Ricci-Tersenghi for useful discussions. 
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